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Abstract. Cosmology contributes a good deal to the investigation of variation of funda- 
mental physical constants. High resolution data is available and allows for detailed analysis 
over cosmological distances and a multitude of methods were developed. The raised demand 
for precision requires a deep understanding of the limiting errors involved. The achievable 
accuracy is under debate and current observing proposals max out the capabilities of todays 
technology. The question for self-consistency in data analysis and effective techniques to 
handle unknown systematic errors is of increasing importance. An analysis based on inde- 
pendent data sets is put forward and alternative approaches for some of the steps involved 
are introduced. 
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1. Introduction 

This work is motivated by numerous findings 
of different groups that partially are in dis- 
agreement witch each other A large part of 
these discrepancies reflect the different meth- 
ods of handling systematic eiTors. Evidently 
systematics are not yet under control or fully 
understood. We try to emphasize the impor- 
tance to take these errors, namely i.e. calibra- 
tion issues, into account and put forward some 
measures adapted to the problem. 
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2. Data 

2.1. Observations 

Q0347-383 is one of the few quasars with 
an absorption spectrum that shows unblended 
strong features of molecular hydrogen. With 
a redshift of z^bs - 3.025 the UV transitions 
are shifted into a spectral range that can be 
observed with earthbound telescopes, render- 
ing it an applicable target for Afi/fj. analysis. 
Hence, Q0347 was su bject of several works 
including the paper by lReinhold et al. I (l2006h 
that ind icated a variation in ix and foHow up pa- 

f-s by King et al. ('2008'), 'Wendt & Reimers] 
)08.) and , Thompson et al.. (.2009.) which re- 
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parabolic fit 
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Fig. 1. The original flux is interpolated by a polyno- 
mial using Neville's algorithm to conserve the local 
flux. 



fitted shift 6.2 mA 
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Fig. 2. Exemplary plot of the sub-pixel cross- 
correlation. The resulting shift is ascertained via 
parabolic fit. In this case the two spectra are in best 
agreement with a relative shift of 6.2 mA. 



evaluated the data and report a result in agree- 
ment with no variation. All recent works on 
Q0347 are based on the same UVES VLT ob- 
serv ationj^ in January 2002 dlvanchik et al. I 
12005 ). For the first time additional observa- 
tional data of Q0347 is taken into account in 
this work. The data were acquired in 2002 at 
the same telescope and recently reduced by 
Paolo MolarcQ. 

2.2. Reduction 

The standard UVES pipeline may not reach 
the desired accuracy for precise determination 
of fundamental physical constants. A major 
source of uncertainty are the calibration spec- 
tra and the limited precision of the available 
lis t of Thorium Argon line s as pointed out 
bv 'Murp hy et al.l(l2008h and [Thompson et all 
(2009). Proper cahbration and data reduction 
will be the key to detailed analysis of potential 
variation of fundamental constants. The influ- 
ence of calibration issues on the data quality 
is hard to measure and the magnitude of the 
resulting systematic error is under discussion. 
This work tries to bypass the estimation of sys- 
tematics and utilizes the uncertainties as part of 
the analysis procedure. 



' Program ID 68.A-0106. 
2 Program ID 68.B-01I5(A). 



2.3. Preprocessing 

The first data set (henceforward referred to as 
set A) consists of nine separate spectra ob- 
served between 7th and 9th of January in 2002. 
The second set of 6 spectra (B) was obtained 
on January 13th in 2002. Prior to further data 
processing the reduced spectra are reviewed in 
detail. 

2.3.1 . Flux error level 

The given error in flux of all 15 spectra was 
tested against the zero level noise in saturated 
areas. A broad region of saturated absorption 
is available near 3913A in the observers frame. 
Statistical analysis revealed a variance corre- 
sponding to ~12Q% of the given error on aver- 
age for the 15 spectra. For each of the spectra 
the calculated correction factor was applied. 

2.3.2. Relative shifts of the 15 spectra 

Due to variation in barycentric velocities and 
grating shifts the individual spectra are sub- 
ject to small shifts - commonly on sub-pixel 
level - in wavelength. To correct for that effect 
all spectra were interpolated by a polynomial 
using Neville's algorithm to conserve the lo- 
cal flux (see Fig. [T]). The resulting pixel step 
on average is 1/20 of the original data. Each 
spectrum was compared to the others. For ev- 
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ery data point in a spectrum the pixel with the 
closest wavelength was taken from a second 
spectrum. Their deviation in flux was divided 
by the quadratic mean of their given errors in 
flux. This procedure was carried out for all 
pixels inside certain selected wavelength inter- 
val^ resulting in a mean deviation of two spec- 
tra. The second spectrum is then shifted against 
the first one in steps of ~ 1 .5 mA. The run of 
the discrepancy of two spectra is of parabolic 
nature with a minimum at the relative shift with 
the best agreement. Fig. |2] shows the result- 
ing curve with a parabolic fit. In this exem- 
plary case the second spectrum shows a shift 
of 6.2 mA in relation to the reference spec- 
trum. The clean parabolic shape verifies the 
approach. All 15 spectra (Ah-B) are shifted to 
their common mean. The average deviation is 
2.3 mA. Section [33] illustrates its influence on 
the data analysis. 

2.3.3. Selection of lines 

The selection of suitable H2 features for the 
final analysis is highly subjective. As a mat- 
ter of course all research groups crosschecked 
their choice of lines for unresolved blends or 
saturation effects. The decision whether a line 
was excluded due to continuum contamina- 
tion or not, however, relied mainly on the ex- 
pert knowledge of the researcher and was only 
partially reconfirmed by the ascertained uncer- 
tainty of the final fitting procedure. This work 
puts forward a more generic approach adapted 
to the fact that we have two distinct observa- 
tions of the same object. 

A selection of 52 lines is fitted separately 
for each dataset of 9 (A) and 6 (B) exposures, 
respectively. In this selection merely blends 
readily identifiable or emerging from equiva- 
lence width analysis are excluded. Each rota- 
tional level is fitted with conjoined line param- 
eters except for redshift naturally. The data are 
not CO added but analyzed simultaneously via 
the fitt ing procedure introduced by Quast et al. I 
(I2005h . For each of the 52 lines there are 



two resulting fitted redshifts or observed wave- 
lengths, respectively, with their error estimates. 
To avoid false confidence, the single lines are 
not judged by their error estimate but by their 
difference in wavelength between the two data 
sets in relation to the combined error estimate. 
Fig.lHshows this dependency. The absolute off- 
set A/leffective to cach Other is expressed in rela- 
tion to their combined error given by the fit: 



^•^effective 



(1) 



■2 +(t] 



Fig. |4] reveals notable discrepancies between 
the two datasets, the disagreement is partially 
exceeding the 5-cr levefl Since the fitting rou- 
tine is known to pro vide pr oper error estimates 
fOuast et al.1 12005') and (IWendt & Reimers I 
2008), the dominating source of error in the 
determination of line positions is due to sys- 
tematic errors. This result indicates calibration 
issues of some significance at this level of pre- 
cision. The comparison of two independent ob- 
servation runs reveals a source of error that can 
not be estimated by the statistical quality of the 
fit alone. For further analysis only lines that 
differ by less than 3 cr are taken into account. 
This criterion is met by 36 lines. Fig. [3] shows 
three H2 features corresponding to the transi- 
tions L5R1, L5P1, L5R2. All have similar sen- 
sitivity towards changes in jd. L5P1 fails the ap- 
plied self consistency check between the two 
data sets and is excluded in the further analy- 
sis. 



3. Results 

For the final analysis the selected 36 lines are 
fitted in all 15 shifted, error-scaled spectra si- 
multaneously. The result of an unweighted lin- 
ear fit corresponds toAyu/// = (15 + 16)x 10"^ 
over the look-back time of ~ 11.5 Gyr for 
Zabs - 3.025. Fig. |5] shows the resulting plot. 



^ Only certain wavelength ranges are taken into 
account to avoid areas heavily influenced by cos- 
mics or areas close to overlapping orders. 



^ Lines fitted with seemingly high precision and 
thus a low error reach higher offsets than lines with 
a larger estimated error at the same discrepancy in 
/iobs- Clearly the lower error estimates merely re- 
flects the statistical quality of the fit, not the true 
value of the specific line position. 
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Fig. 3. Part of the co added spectrum near 4176A. 
The data however, were not co added for the fit. 
L5R1 and L5R2 match the 3-cr criterion, L5P1 does 
not. 
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Fig. 5. The unweighted fit for Q0347 corresponds 
to Ap/yu = (15 + 16) X 10"^. Note, that at such a high 
scatter ZKi=o differs from z by less than 1 cr^. 
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Fig. 4. Selection of 52 seemingly reasonable lines 
to be fitted separately for each dataset of 9 and 6 ex- 
posures, respectively. Their absolute offset A/ldfcctive 
to each other is expressed in relation to their com- 
bined error given by the fit (see Eq[T](. 



Note the comparable large scatter in deter- 
mined redshift when combining the two in- 
dependent observations. The approach to ap- 
ply an unweighted fit is a consequence of the 
unknown nature of the prominent systematics. 
Uncertainties in wavelength calibration can not 
be expressed directly as an individual error per 
line. For this work the calibration errors and 
the influence of unresolved blends are assumed 
to be dominant in comparison to individual fit- 
ting uncertainties per feature. For the following 
analysis the same error is adopted for each line. 
With an uncertainty in redshift of 1 x 10"^ we 
obtain: Aju/;U = (15 ± 6) x 10"^. 



However the goodness-of-fit is below 1 ppm 
and is not self consistent. Judging by that and 
Fig. |5l a reasonable error in observed redshift 
should at least be in the order of 5 x 10"^. The 
weighted fit gives: A;u/// - (15 + 14) x 10"*. 



3. Uncertainties in the sensitivity 
coefficients 

At the current level of precision, the influ- 
ence of uncertainties in the sensitivity coef- 
ficients Ki is minimal. It will be of increas- 
ing importance though when wavelength cal- 
ibration can be improved by pedantic demands 
on future observations and eventually Laser 
Frequency Comb calibration will allow for 
practically arbitrary precision and uncertain- 
ties in the calculations of sensitivities will play 
a role. Commonly the weighted fits neglects 
the error in Ki. The merit function for the 
case of a straight-line fit with errors in both co- 
ordinates is given by: 

X {a,b)^ ) — ^ (2) 

1=0 .V/ XI 

and can be so lved numerically with valid ap- 
proximations dLybanon II 19841) . At the current 
level even an error in K, of about 10% has 
merely impact on the error estimate in the or- 
der of 10 ^ 

Alternatively the uncertainties in Ki can be 
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translated into an uncertainty in redshift via 
the previously fitted slope. The results of this 
ansatz are similar to the fit with errors in both 
coordinates and in general this is simpler to im- 
plement. 

Another possibility is to apply a gaussian error 
to each sensitivity coeflicient and redo the nor- 
mal fit a few dozen times with alternating vari- 
ations in Ki. Again, the influence on the error- 
estimate is in the order of 1 ppm. 
The different approaches to the fit allow to es- 
timate its overall robustness as well. 

3.2. Individual line pairs 

A/u/fi can also be obtained by using merely 
two lines that show different sensitivity to- 
wards changes in the proton-to-electron mass 
ratio. Another criterion is their separation in 
the wavelength frame to avoid pairs of lines 
from different ends of the spectrum. Several 
tests showed that a separation of AA < 11 OA 
and a range of sensitivity coefficients K1-K2 > 
0.02 produces stable results that do not change 
any further with more stringent criteria. Note, 
that pairs spread across two neighboring orders 
(~ 50A) show no striking deviations. Fig]?] 
shows the different values for Afi/fi derived 
from 52 line pairs that match the aforemen- 
tioned criteria. Note, that a single observed line 
contributes to multiple pairs. The gain in sta- 
tistical signific ance by this sorting i s limited as 
pointed out by lMolaro et al. I (l2008h . The large 
scatter favors the median as result which pro- 
duces Afi/fu = 13 X 10"^. The scatter is then 
related to uncertainties in the wavelength de- 
termination which is mostly due to calibration 
errors. The standard error is 8 x 10"^. 

3.3. Influence of the preprocessing 

Section 12.3.21 describes the initial shift to a 
common mean of all 15 spectra. The com- 
plete analysis was redone with error-scaled but 
unshifted spectra and the ascertained line po- 
sitions of both runs compared. Fig. |6] shows 
the difference for each H2 line in mA over the 
corresponding sensitivity coefficients Kj. The 
plotted line is a straight fit. Cleai^ly the slope is 
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Fig. 6. Variation in fitted positions for all lines with 
and without initial correction for shifts in between 
the 15 spectra. The slope of the fit is dominated by 
three lines. 
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Fig. 7. Afi/jA derived from individual line pairs (52) 
which are separated by less than 11 OA and show a 
difference in sensitivity of more then 0.02. 
The median corresponds to A/i/fi = 13 x 10"*. 



dominated by three individual lines whose fit- 
ted centroids shifted up to 5.5 mA due to the 
preprocessing. These three lines in particular 
produce a trend towards variation in fi when 
grating shifts and other effects are not taken 
into account. This single-sided trend probably 
occurred by mere chance but at such low statis- 
tics it influences the final result. 
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